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We report systematic Cu- and F-NMR measurements of five-layered high-T c cuprates 
Ba2Ca4Cu50io(F,0)2- It is revealed that antiferromagnetism (AFM) uniformly coexists with su- 
perconductivity (SC) in underdoped regions, and that the critical hole density p c for AFM is ~ 
0.11 in the five-layered compound. We present the layer-number dependence of AFM and SC phase 
diagrams in hole-doped cuprates, where p c for n-layered compounds, p c (n), increases from p c (l) 
~ 0.02 in LSCO or p c (2) ~ 0.05 in YBCO to p c (5) ~ 0.11. The variation of p c (n) is attributed 
to interlayer magnetic coupling, which becomes stronger with increasing n. In addition, we focus 
on the ground-state phase diagram of Cu02 planes, where AFM metallic states in slightly doped 
Mott insulators change into the uniformly mixed phase of AFM and SC and into simple d-wave 
SC states. The maximum T c exists just outside the quantum critical hole density, at which AFM 
moments on a Cu02 plane collapse at the ground state, indicating an intimate relationship between 
AFM and SC. These characteristics of the ground state are accounted for by the Mott physics based 
on the t-J model; the attractive interaction of high-T c SC, which raises T c as high as 160 K, is an 
in-plane superexchange interaction Ji n (~ 0.12 eV), and the large Ji„ binds electrons of opposite 
spins between neighboring sites. It is the Coulomb repulsive interaction U (> 6 eV) between Cu-3d 
electrons that plays a central role in the physics behind high-T c phenomena. 



I. INTRODUCTION 

Despite extensive research for more than a quarter 
of a century, there is still no universally accepted the- 
ory about the mechanism of superconductivity (SC) in 
copper-oxide superconductors. The main controversy ex- 
ists over the attractive force that forms Cooper pairs, 
which leads to a remarkable high SC transition temper- 
ature T c . It has been known that the hybridization be- 
tween Cu(3d x 2_ y 2) and 0(2pa) orbitals induces a large 
in-plane superexchange interaction Ji n ~ 0.12 eV or 1300 
K in nondoped Cu0 2 planes, and that strong antifer- 
romagnetic (AFM) correlations exist even in SC states. 
Experiments have observed local magnetic moments in 
metallic or SC regions [iT-fl3j. and theories have pointed 
out a strong relationship between AFM and SC in under- 
doped regions [l4T-l3lT|. Those reports have shown that 
AFM plays a key role in the mechanism of high-T c SC in 
cuprates. 

We have reported on AFM and SC properties in mul- 
tilayered cuprates [32h35| | , where long-range static AFM 
orders uniformly coexist with SC in underdoped regions 
[3~jj . On the other hand, the uniform coexistence has not 
been observed in typical high-T c cuprates such as single- 
layered La2_xSr x CuC"4 (LSCO) [h| HH and bi-layered 
YBa 2 Cu 3 6+y (YBCO) @, G3. The difference of the 
phase diagrams in various cuprates is possibly attributed 
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to the strength of interlayer magnetic coupling, which 
would become stronger with an increase in the stacking 
number n of Cu02 layers in a unit cell. In addition to 
that, it is expected that disorder effects, which is in asso- 
ciation with chemical substitutions, are enhanced when 
n is small pjj EH H3, HI] , and that such a disorder masks 
the intrinsic nature of Cu02 planes. In order to under- 
stand the relationship between AFM and SC, therefore, it 
is necessary to investigate the n dependence of electronic 
properties in underdoped regions. 

Apical-fluorine multilayered high-T c cuprates 
Ba 2 Ca„_iCu„0 2 „(F,0) 2 (02(n-l)nF) provide us 
with the opportunity to investigate the n-dependent 
AFM-SC phase diagrams. 02(n-l)nF comprises a stack 
of n Cu0 2 layers, as shown in Figs. IUa)-njd), and 
consists of inequivalent Cu0 2 layers: an outer Cu0 2 
plane (OP) in a five-fold pyramidal coordination and an 
inner Cu0 2 plane (IP) in a four-fold square coordina- 
tion. The substitution of oxygen 2 ~ for apical fluorine 
F 1_ results in doping hole carriers, increa sing T r from 
underdoped to optimally-doped regions [33J, [35|, |39|-|4l| . 

In this paper, we report the AFM and SC phase 
diagram in five-layered Ba 2 Ca4Cu 5 Oi (F,O) 2 (0245F) 
by means of 63 Cu- and 19 F-NMR, and present the n- 
dependence of the phase diagram from n=l to 5. We 
highlight the fact that the ground-state phase diagram 
derived from the present study is in good agreement with 
theoretical predictions based on the t-J model. This re- 
sult supports the idea that the in-plane superexchange in- 
teraction J in plays a vital role as the glue to form Cooper 
pairs or mobile spin-singlet pairs. 

The contents of this paper are as follows: after ex- 
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FIG. 1. (color online) Crystal structures of Ba 2 Ca„_iCu„0 2 n(F,0)2 with (a) n=2 (0212F), (b) n=3 (0223F), (c) n=4 (0234F), and 
(d) n=5 (0245F). Here, n is the number of stacked CUO2 layers. The heterovalent substitution of F 1_ for O 2- at apical sites decreases 
the hole density p. Outer and inner Cu02 planes are denoted as OP and IP, respectively. Note that, in a precise sense, there are two 
kinds of IP layers: the innermost inner CUO2 plane (IPq) and the two adjacent inner ones (IPi). 



perimental details in Section II, we provide experimental 
results and discussions; in Section III, we show system- 
atic Cu- and F-NMR measurements on four 0245F sam- 
ples, which is evidence of the uniform coexistence of AFM 
and SC; in Section IV, we present the n-dependence of 
the AFM and SC phase diagrams in 02(n-l)nF, and fi- 
nally construct the ground-state phase diagram inherent 
in hole-doped cuprates. 



II. EXPERIMENTAL DETAILS 

Polycrystallinc powder samples used in this study were 
prepared by high-pressure synthesis, as described else- 
where [39|, • Powder X-ray diffraction analysis shows 
that the samples comprise almost a single phase. As dis- 
cussed later, the sharp Cu-NMR spectral width at T=300 
K assures the good quality of the samples. We prepared 
four samples of Ba2Ca4Cu50io(Fj / Oi_ 2/ )2 (0245F) and 
determined the values of T c by the onset of SC diamag- 
netism using a dc SQUID magnetometer. The four 0245F 
samples exhibit a systematic decrease in T c , as listed in 
Table HI as the nominal amount of fluorine F 1_ (i.e., y) 
increases. The heterovalent substitution of F 1_ for O 2- 
at apical sites (see Fig. [T]) decreases the hole density p in 
apical-F compounds [33l l35l. [39l - l4lT | . Note that it is dif- 
ficult to investigate the actual fraction of F 1_ and O 2- 
[13, 1!!!. For NMR measurements, the powder samples 
were aligned along the c axis in an external field H ex of 
~ 16 T and fixed using stycast 1266 epoxy. NMR ex- 
periments were performed by a conventional spin-echo 
method in the temperature (T) range of 1.5 — 300 K. 



III. RESULTS 

A. Cu-NMR 

1. 63 Cu-NMR spectra and estimation of nuclear quadrupole 
frequency 63 vq 

Figures (D(a), gjb), gfc), and^d) show typical 63 Cu- 
NMR spectra of the central transition (1/2 <=> —1/2) 
for 0245F((J1), 0245F(}J2), 0245F(p), and 0245F(jj4), re- 
spectively. The field-swept NMR spectra were measured 
with H ex perpendicular to the c axis {H ex _!_ c), and the 
NMR frequency o>o was fixed at 174.2 MHz. As shown 
in Fig. [ljd), 0245F has two kinds of CuC>2 planes: an 
outer CuC>2 plane (OP) and an inner plane (IP). There- 
fore, the two peaks in the NMR spectra correspond to 
OP and IP. The assignment of NMR spectra to OP and 
IP has been already reported in previous literature [44r - 
|46j. Here, note that Cu-NMR spectra for the innermost 
IP (IP in Fig. Old)) and the two adjacent IP (IPi in 
Fig. Did)) overlap each other, which suggests that their 
local doping levels are not so much different [4l| . Hence- 
forth, we do not distinguish IPo and IPi in this paper. In 
0245F(jjl), the Cu-NMR spectra for both OP and IP are 
observed at T=280 K, whereas the IP's spectrum disap- 
pears at low temperatures, as shown in Fig. Ufa). This is 
because AFM correlations develop upon cooling as in the 
case of three-layered 0223F (35[ and four-layered 0234F 
[33 | . We have also reported on the loss of Cu-NMR in- 
tensity due to spin dynamics in a previous paper [32j |. 
In 0245F(jj2), the Cu-NMR spectra are observed for OP 
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TABLE I. List of physical properties in Ba2Ca4Cu50io(F !/ Oi_ a )2 (0245F) used in this study. The T c values were determined 
by the onset of SC diamagnetism using a dc SQUID magnetometer. Here, note that the amount of F 1 ", y, is nominal. It is 
difficult to determine the actual fraction of F _1 and O 2- at the apical sites [iQ, 43;]. The hole density p, the nuclear quadrupole 
frequency 63 uq, and the full- width at half-maximum (FWHM) are separately evaluated for OP and IP from the Cu-NMR 
measurements at T=300 K (see text). 
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FIG. 2. (color online) Typical 63 Cu-NMR spectra for (a) 
0245F()tl), (b) 0245F(P), (c) 0245F(P), and (d) 0245F((|4). The 
spectra were measured with H ex _L c and at wo=174.2 MHz. For 
(a) 0245F((tl) and (b) 0245F(t)2), the spectra of IP disappear at low 
temperatures due to the development of antiferromagnetic correla- 
tions. For (c) 0245F())3) and (d) 0245F(t(4), no signal is obtained 
for IP even at T=300 K. 



and IP at T=300 K, as shown in Fig. H». When de- 
creasing T, however, not only the spectrum of IP but also 
that of OP disappears at low temperatures. Moreover, 
in 0245F(U3) and 0245F(|)4), the IP's spectra are not ob- 
served even at T=300 K. These marked differences in 
the NMR spectra among the four samples suggest that 
AFM correlations become stronger as p decreases from 
0245F()J1) to 0245F()J4). The values of p for the present 
samples, listed in Table HI are discussed later. Here, note 
that the full- width at the half maximum (FWHM) of the 
NMR spectra at T=300 K decreases from 0245F(jjl) to 
0245F(jj4), as presented in Table fl] This is because the 
disorder associated with the atomic substitution is re- 
duced with increasing y, i.e., decreasing the amount of 
the O -2 substitution at the apical-F sites. In the same 
sample, the NMR spectral width for OP is much broader 



than that for IP; OP is closer to the Ba-F layer (see Fig. 
[1]), which is the source of the disorder due to the atomic 
substitution at apical-F sites. The values of FWHM for 
OP are less than ~ 300 Oe, which is narrower than those 
for Bi- and Tl-compounds [13, |48| . On the other hand, 
FWHM for IP is less than ^110 Oe, which is comparable 
to or even narrower than those for Y- and Hg-compounds 
[4^ | . These NMR linewidths point to the good quality of 
the present samples. 



0.8 



0.6 



0.4 



1 1 1 
• (#1)OP 


1 < 1 


- a (#2) OP 




■ (#3) OP 




+ (#4) OP 


(#1)IP 




" 7 = 300 K 

1 . 1 


A (#2) IP 



4 6 
lYHres)' 2 (10" 5 MHz" : 



FIG. 3. (color online) NMR frequency a;o-dependence of H reB 
plotted as (ui - 7jv-H r es)/7iv#res vs (7j\r#res) -2 . According to 
Eq. |[T}, the values of 63 vq are estimated as listed in Table U 

According to the second-order perturbation theory for 
the nuclear Hamiltonian with H ex ± c [5(1 [H| , the NMR 
shifts of the spectra in Fig. [5]consist of the Knight shift K 
and the second-order quadrupole shift. The NMR shifts 
are expressed as 



w — lNH re 



= K 



1 



(1) 



l N H res " ' 16(l + K)(j N H res ) 2 

where 7^ is a nuclear gyromagnetic ratio, H res is an 
NMR resonance field, and vq is a nuclear quadrupole 
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frequency. In order to estimate 63 for the present 
samples, we have measured the uj dependence of H res 
at T=300 K in a range of 110.2 - 176.2 MHz. The 
obtained data set of luq and H res is plotted as (ujq — 
"f N H res )/j N H res vs (jNHresy 2 in Fig. OS Based on 
Eq. ([1]), we estimated 63 vq from the slope of the linear 
line in the figure. The obtained 63 vq values, listed in 
Table HI are comparable with those in other multilayered 
cuprates [U 0, H| El HI ■ As shown in Table U %z v Q 
decreases from 0245F(j}l) to 0245F(|j4); this reduction of 
63 v Q shows that p decreases from 0245F(j}lVto 0245F(t)4), 
as observed in other hole-doped cuprates fill |5~4T - [58T | . 



Cu-NMR shift and estimation of hole density p 
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FIG. 4. (color online) T-dependences of 63 Cu Knight shift K S (T) 
with H ex i_ c for (a) 0223F (cited from Ref. gj), (b) 0245F(J1), 
(c) 0245F(tt2), (d) 0245F((t3), and (e) 0245F(|)4). 

According to Eq. ((T|), the Knight shift K for the Cu- 
NMR spectra in Fig. [2] is estimated by subtracting the 
second order quadrupole shift from the total NMR shift. 
In order to estimate K, we use the Q3 vq values listed 
in Table H Here, 63 vq(IP) ~ 7.0 to 7.5 MHz is assumed 
for 0245F(p) and 0245F(jj4) since 63 v Q becomes smaller 



with decreasing p. In high-T c cuprates, K comprises a In- 
dependent spin part K S {T) and a T- independent orbital 
part K or b as follows: 



K = K S (T)+K 



(2) 



The T-dependences of K S (T) with H ex _L c for 0245F(jjl), 
0245F(j}2), 0245F(j}3), and 0245F(|4) are displayed in 
Figs. HJb), 0Jc), HJd), andHJe), respectively. Here, K orb 
was determined as ~ 0.22 %, assuming K S (T) ~ at a 
T = limit. 

As shown in the figures, the room temperature value 
of K S (T) decreases with decreasing p from 0245F(jJl) to 
0245F()J4). The values of p(OP) and p(IP), which are 
summarized in Table U are separately evaluated by us- 
ing the relationship between K s (300 K) and p [41j |. The 
quantity p(lP) is smaller than p(OP) because IP is far 
from charge reservoir layers, which has been usually ob- 
served in multilayered cuprates 0, |4(| [5l|. As for IP 
in 0245F(jl3) and 0245F(JJ4), it is impossible to directly 
estimate p(lP) from ^(300 K). Therefore, we tenta- 
tively estimate p(lP) from Fig. [5j Figure [5] shows p(OP) 
and p(IP) as functions of the average hole density p av 
for 0245F and another five-layered compound Hgl245 
[12. Here, p av is defined as p av = {2p(OP)+3p(IP))/5. 
As shown in the figure, p(OP) and p(lP) systemati- 
cally decrease with the reduction of p av , which allows 
us to extrapolate p av for 0245F(jj3) and 0245F()}4) from 
p(OP). As for 0245F(j}3), p atI =0.065 is expected from 
p(OP)=0.083 on the linear line for p(OP) vs p av . Fur- 
thermore, the plot for p(lP) vs p av gives a tentative value 
of p(IP)=0.053 from p at ,=0.065. The value of p(IP) for 
0245F(t)4) is also obtained by adopting the same proce- 
dure. As summarized in Table HI p(OP) and p(lP) de- 
crease with the increase of the fluorine content ?/, which 
reduces T c systematically. 

The T-dependence of K s (T) for the four 0245F samples 
is different from that for paramagnetic superconductors, 
suggesting a possible AFM order at low temperatures. 
As an example of a paramagnetic multilayered com- 
pound, we show in Fig.rjja) the T-dependence of K S (T) 
for optimally-doped three-layered Ba2Ca2Cu30g(F,0)2 
(0223F) with T c =120 K [HED]. As shown in Fig. [Da), 
K S (T) decreases upon cooling down to T c in association 
with the opening of pseudogap HU . The steep de- 
crease of K s (T) below T c is evidence of the reduction in 
spin susceptibility due to the formation of spin-singlet 
Cooper pairing. These behaviors in K S (T) are common 
in hole-doped cuprates that are paramagnetic supercon- 
ductors. On the other hand, in the case of 0245F(fll), 
K s (T) in Fig. 0Jb) shows a T-dependence totally differ- 
ent from that in Fig. |4|a); K S (T) for IP can not be de- 
termined below T ~ 200 K because of the disappearance 
of Cu-NMR spectra, and K s (T) for OP shows an upturn 
at T - 85 to 90 K. These unusual behaviors of K 8 {T) 
suggest a possible AFM transition at IP in 0245F(jJl) 
with the Neel temperature Tv ~ 85 K. A similar T- 
dependence in K S (T) has been reported in optimally- 
doped five-layered Hgl245 jHIH, where K S (T) for OP 
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FIG. 5. (color online) Plot of p(OP) and p(IP) as func- 
tion of p av for 0245F(ftl), 0245F(Jt2), 0245F(fl3), and 0245F()|4). 
Here, p a v is the average value of p(OP) and p(IP), denned as 
p a „ = (2p(OP)+3p(IP))/5. The data for Hgl245 is cited from 
Ref. gg. As for IP in 0245F(t(3) and (jj4), it is impossible to di- 
rectly estimate p(IP) from K s (300 K). Therefore, we tentatively 
estimate p(IP) from the linear lines in the figure. As shown in 
the figure, p(OP) and p(IP) systematically decrease with the re- 
duction of pan, which allows us to extrapolate p av for 0245F(()3) 
and 0245F(t(4) from p(OP). As for 0245F(B3), p a „=0.065 is ex- 
pected from p(OP)=0.083 on the linear line for p(OP) vs p a v 
Furthermore, the plot for p(IP) vs p a v gives a tentative value of 
p(IP)=0.053 from p a „=0.065 as shown by an open circle. The 
value of p(IP) for 0245F(J4) is also obtained by adopting the same 
procedure. As summarized in Table I, p(OP) and p(IP) decrease 
with the increase of the fluorine content y, which reduces T c sys- 
tematically. 



shows an upturn at T ~ 55 K with the AFM transition 
at IP. In 0245F(jj2), the values of K S (T) are unmeasur- 
able for both OP and IP at low temperatures as shown 
in Fig.glc). In 0245F(jJ3) and 0245F(J4), K S (T) for OP 
determined only at high temperatures, and K S (T) for IP 
is not obtained in all measured T-ranges, as shown in 
Figs. Bid) andBJe). These unusual behaviors of K S (T) 
suggest that AFM orders occur at both OP and IP in 
0245F(jj2), 0245F(p) and 0245F(jj4), and that T N in- 
creases with decreasing p from 0245F(JJ2) to 0245F(j}4). 



B. Evidence for AFM order probed by Cu-NQR 
and zero-field Cu-NMR 

As shown in Fig. H the 63 Cu-NMR spectra for 0245F 
are lost when T decreases. This is due to the marked 
development of AFM correlations, which suggests that 
static AFM orders occur at low temperatures. NMR 
measurements at H ex —0 T sensitively detect evidence of 
static AFM orders, as explained below. 

In general, the Hamiltonian for Cu nuclear spins (/ = 
3/2) with an axial symmetry is described in terms of the 
Zeeman interaction T-Cz due to a magnetic field i/, and 
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FIG. 6. (color online) (a) Cu-NQR spectrum at T=1.5 K for 
three-layered 0223F with T c =120 K (cited from Ref. ;34]). The 
spectrum is typical of paramagnetic multilayered compounds, (b)- 
(e) Cu-NQR or zero-field NMR spectra at T=1.5 K for 0245F. The 
data for 0245F(J)4) in (e) is cited from Ref. 13411 . The bars represent 
resonance frequencies estimated by using Eq. (|3jl with vq values in 
Tableland H ini values in Table[ll] The curves in the figure are Cu- 
NMR spectral simulations for OP and IP based on the positions 
of the bars, which represents resonance frequencies estimated by 
using Eq. and simulations to the total Cu-NMR spectra. 



the nuclear-quadrupole interaction T~Lq as follows: 
n = H z + H Q 

e 2 qQ 



- lN hi-H- 



41(21 - 1) 



(3/^,-/(7+1)), (3) 



where eQ is the nuclear quadrupole moment, and eq is 
the electric field gradient at a Cu nuclear site. In "Hq, 
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the nuclear quadrupole resonance (NQR) frequency is de- 
fined as vq — e 2 qQ/2h. In paramagnetic substances, an 
NQR spectrum is observed due to the second term in 
Eq. §3§ when H=H ex —0 T. On the other hand, in mag- 
netically ordered substances, an internal magnetic field 
Hi n t is induced at Cu sites; in addition to the second 
term, the first term in Eq. contributes to the nu- 
clear Hamiltonian even if H ex —0 T. Therefore, the onset 
of a magnetically ordered state is observed as a distinct 
change of the spectral shape at H ex =0 T. 

Figure |6ja) shows the Cu-NQR spectrum of three- 
layered 0223F with T c =120 K [34[ as an example of para- 
magnetic multilayered cuprates. The two peaks corre- 
spond to OP and IP, and both peaks include two com- 
ponents, 63 Cu and 65 Cu. It is assured that 0223F is a 
paramagnetic superconductor because an NQR spectrum 
is obtained for OP and IP at H ex ~0 T. Here, note that 
an NQR spectrum similar to that in Fig. H^a) should be 
obtained at H ex —0 T when a measured material is para- 
magnetic. The spectra obtained at H ex =0 T for five- 
layered 0245F are, however, totally different from that of 
paramagnetic 0223F as shown in Figs. EfbJ-Ete). 

Figure [B^b) shows the spectrum measured at H ex =0 T 
and T=1.5 K for 0245F(jjl). An NQR spectrum is ob- 
served for OP, and the NQR frequency 63 vq ~ 14 MHz 
corresponds approximately to the value estimated from 
Fig. [3] (see Tabic [I]). On the other hand, the spectrum for 
IP is different from the NQR spectrum for IP in 0223F; 
the resonance frequency of IP in 0245F((|1) is significantly 
larger than 63 vq ~ 8.3 MHz, which was estimated from 
Fig. El According to Eq. @, resonance frequencies in- 
crease when Hi n t is induced at Cu sites by the onset of 
AFM orders. The bars in Fig. [BJb) represent the reso- 
nance frequencies estimated by using Eq. (O on the as- 
sumption of H^t = T for OP and of H mt ~ 2.3 T 
(_L c) for IP. This reveals that Hi nt ~ 2.3 T is induced 
at IP by spontaneous AFM moments Mafm due to the 
AFM order. As shown in the figure, there are 63 Cu and 
65 Cu components, and each Cu component has one center 
peak and two satellite peaks when Hint ^ 0. However, 
due to the poor frequency resolution that is related to 
a weak signal to noise ratio, those signals are not well 
resolved. The weak resolution may be also attributed to 
the inhomogeneities of vq and of the size or direction of 
Mafm- Here, we tentatively represent spectra of the IP 
by a single Lorentzian as shown in Fig. [5] (b). 

Figures|6jc),[6jd), and[6je) show the spectra measured 
at H ex =0 T and T=1.5 K for 0245F(jJ2), 0245F(j}3), and 
0245F(fl4), respectively. In the three samples, the spec- 
tra for both OP and IP are totally different from the 
NQR spectra for 0223F shown in Fig. EJa), which sug- 
gests AFM orders at both OP and IP. With decreasing 
p from 0245F(jJ2) to 0245F(j}4), the spectrum shifts to 
higher frequency regions, as shown in the figures. As in 
the case of IP in 0245F(jJl), the bars in the figures repre- 
sent the resonance frequencies estimated by using Eq. ([3]) . 
Here, we used the values of 63 vq listed in Table Q] and as- 
sumed the values of Hi„ t listed in Table |TTJ The quantity 



H^t is converted into Mafm, as listed in Table [TJl by us- 
ing the relation of H int = \A hf \M A FM = \A-AB\M A fm- 
Here, A and B are the on-site hyperfine field and the su- 
pertransferred hyperfine field, respectively. A ~ 3.7 and 
B ~ 6.1 T/fJ>B are assumed, which are typical values in 
multilayered cuprates in underdoped regions [ilf . Here, 
we comment on some distributions of each resonance fre- 
quency presented in Fig. |BJ As mentioned above, the 
weak resolution may be due to the inhomogeneities of uq 
and Hint- However, a possibility of spin-glass states is 
ruled out in these compounds. It is reported that the 
spin-glass phases cause Hint to largely distribute in as- 
sociation with random directions of frozen magnetic mo- 
ments, which does not allow us to obtain zero-field spec- 
tra in limited frequency regions. In fact, 19 F-NMR stud- 
ies in the following section assure onsets of long-range 
three dimensional AFM orders. 

TABLE II. List of H int , Mafm, and T N for 0245F. We used 
Hint values presented below in order to obtain the bars in 
Fig. |6j The values of Mafm are deduced by using the re- 
lation of Hint = \A hf \MAFM = \A — 4B\Mafm, where A 
and B are the on-site hyperfine field and the supertransferred 
hyperfine field, respectively. A ~ 3.7 and B ~ 6.1 T//j,b are 
assumed, which are ty pica l values in multilayered cuprates 
in underdoped regions [44J. The values of Tjv are estimated 
from the upturn in K S (T) for 0245F(|J1) and from F-NMR 
measurements for 0245F(|}2), 0245F(jj3), and 0245F(jj4) (see 
text) . 
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C. F-NMR 

1. 19 F-NMR spectra 

It is difficult to deduce the T dependences of 
Mafm (OP) and Mafm (IP) from the zero-field Cu-NMR 
measurements because the nuclear spin relaxation at Cu 
sites is enhanced greatly as T approaches Tjv. Instead, 
19 F-NMR is used to probe the T dependence of the in- 
ternal field Hi ri t(F) at apical-F sites, which is induced 
by Mafm- Figure [7] shows the T dependences of 19 F- 
NMR spectra obtained by sweeping frequencies with H ex 
(=4.245 T) parallel to the c-axis. 

Figure [Jd) shows the T dependence of the F-NMR 
spectrum for 0245F((j4), which has the lowest T c of the 
present samples. A sharp spectrum is observed with 
a single peak at T = 240 K, but the spectrum splits 



7 




FIG. 7. (color online) T dependences of 19 F-NMR spectra for (a) 0245F(J1), (b) 0245(^2), (c) 0245F((|3), and (d) 0245F(j|4). Here, H sx 
is parallel to the c axis and is fixed at 4.245 T. The spectra in (b)-(d) split into two at low temperatures, pointing to the onset of AFM 
orders. Note that the spectrum for (a) 0245F((J1) exhibits a single peak in all measured T ranges because OP is paramagnetic (see text). 
The dashed line indicates 19 if c =0, and the solid lines are spectral simulations to estimate ui, \Hi n t e(F)|i an d 1S K C . 



into two peaks at low temperatures. As discussed be- 
fore, 0245F(j}4) shows the AFM order with M A fm (OP) 
- 0.14 and M AFM (IP) ~ 0.20 at T=1.5 K. Therefore, the 
spectral splitting suggests the development of Hi nt (F) in- 
duced by Mafm- Figure EJc) presents the T dependence 
of the resonance frequency ui in the 19 F-NMR spectra, 
which is deduced through spectral simulations shown by 
solid lines in Fig. [7£d). As shown in Fig. [5Jc), the F- 
NMR spectra split below T=175 K, suggesting the AFM 
ordering below the Neel temperature K. 

Figures [7fb) and [Tfc) show the T-dependences of the 
F-NMR spectra for 0245F(J}2) and 0245F(j}3), respec- 
tively. As in the case of 0245F(j}4), the F-NMR spectra 
split into two at low temperatures. Below T c =75 K in 
0245F()j2) and T c =65 K in 0245F(j}3), however, the two 
spectral peaks are blurred due to the spectral broaden- 
ing related to vortex states in the SC phase. Note that 
M AFM in 0245F(j}2) and 0245F(jj3) is smaller than that 
in 0245F(jJ4), and that the vortex- related spectral broad- 
ening in 0245F(t)2) and 0245F(j}3) is larger than that in 
0245F(t)4) due to the higher T c values. Therefore, in con- 
trast to 0245F()J4), it is difficult to separate the two peaks 
in 0245()J2) and 0245F(j}3) below T c . The T-dependences 
of w for 0245F(jj2) and 0245F(jj3) are shown in Figs.^a) 
and [5th), suggesting the AFM transitions at T/v ~ 95 K 
and ~ 145 K, respectively. 



Figure [TJa) shows the T-dependence of the F-NMR 
spectrum in 0245F(JJ1), which has the highest T c in the 
present samples. The spectral shape is unchanged upon 
cooling above T c =85 K; the spectrum shifts to lower fre- 
quency regions due to the SC transition below T c . As 
shown in Fig.^b) and listed in TableHH OP in 0245F()J1) 
is paramagnetic at T=1.5 K, whereas IP shows an AFM 
order with Mafm@P) ~ 0.11 \ib- Therefore, the fact 
that the spectrum does not split at any temperature 
shows that the F-NMR spectra with H ex parallel to the 
c axis is not affected by MafmJJP), as we already dis- 
cussed in the previous report [341 ] . 

2. T dependence of internal field at apical-F site 

The resonance frequency uj of 19 F-NMR with H ex par- 
allel to the c-axis is expressed by 

U> ~ W lN H ex (1 + 19 K C ) ± 19 7 Ar|^nt,c(F)|, (4) 

where Hi nt>c (F) is the component of Hi nt (F) along the 
c-axis, 19 7at is the 19 F nuclear gyromagnetic ratio, and 
19 K C is the Knight shift. The plus (minus) sign of 
\Hint,c(F) \ in Eq. (0| corresponds to its parallel (antipar- 
allel) component along the c-axis. As shown in Figs.[5Ja)- 
|S£c), the spectral splitting Aw increases below Tjv due 



8 



N 



O 



1 70.08 
1 70.04 

170.00 

8 
6 
4 
2 


0.02 
0.00 
-0.02 





(a)0245F(#2) 

5§x> o o o o o o o. 






( 


(d) 0245F(#2) 
T N =95 K 

o 

✓">. ....... 


i 

o 
o 
o 
o 

rP. . 


rt£D O O O O O O O 
f 

Tc=75 K 

(g)0245F(#2) 



(b) 0245F(#3) 



o o o 



(e) 0245F(#3) 
r N =145 K 





(c) 0245F(#4) 



O O o o 



°o ! 


(f) 0245F(#4) 






7n=175K 



O O o o 



Tc=52 K 



(i) 0245F(#4) 



50 100 150 200 250 300 50 100 150 200 250 300 50 100 150 200 250 300 

T(K) r(K) T(K) 



FIG. 8. (color online) (a)-(c) T dependence of a; for 0245F(J2), 0245F(ft3), and 0245F(J4). The values of u are estimated from spectral 
simulations in Fig. [7J For (a) 0245F(jJ2) and (b) 0245F(()3), it is impossible to determine ui below T c owing to the marked spectral 
broadening, which is in association with the distribution of vortices in the SC mixed state, (d)-(f) T dependence of |//i n t,c(F)| estimated 
from Aaj=2x 19 7jv|-Hint,c(F)|. Here, Aoj is the width of the spectral splitting in F-NMR spectra. The solid line is -Hi n t,c(F)oc Mafm(T) = 
M AFM (0)(1 - T/T N ) - 5 , which is in good agreement with the experiment down to T c . (g)-(i) T dependence of Knight shift 19 K C . The 
decrease in 19 K C below T c is associated with SC diamagnetic shifts. The data in (c), (f), and (i) for 0245F(fl4) are cited from Ref. [34ll . 



to the development of flj„t )C (F), which is in proportion 
to Mafm(OP)- If the direction of Mafm (OP) were in 
the basal Cu02 planes with a wave vector Q = (77, 7r), 
Hint,c(F) would be cancelled out at the apical-F site; 
the F ions are located at a magnetically symmetric posi- 
tion. We have reported that H int ^{¥) is induced because 
Mafm {OP) is directed out of the planes within a few de- 
grees of the canting angle [34| . Refer to Ref. [34[ for more 
detailed discussions on the origin of |if,* ni)C (F)|. 

Figure [S^f) shows the T-dependence of 
\H int , c (F)\ for 0245F(jj4), which is obtained from 
Aw=2x 19 7jv|.ffi TltjC (F)|. In order to focus on the 
T evolution of Mafm, Mafm{T), the relationship 
H int , c {F) oc Mafm (T) = M AFM (0)(1 - T/T N )°- 5 is 
displayed as the solid line in Fig. |UJf). This power-law 
variation of Mafm(T) is in good agreement with the 
experiment down to T c — 52 K, which suggests the 
three-dimensional long-range order of Mafm- The 
same T-dependence of Mafm(T) has been reported in 
slightly-doped LSCO compounds which exhibit AFM 
ground states [62| . It is also possible to reproduce the T 
dependence down to T c =52 K by another relationship 
H mt , c (F) cx Mafm (T) = M AFM {0){1 - (T/Tjv) 3 / 2 ) - 5 
(not shown), which has been theoretically predicted 
on an itinerant AFM metal [63[- In any case, the T 
evolution of Mafm assures that three dimensional 
AFM orders set in below Tjy. Figure HJe) shows the 
T-dependence of |i2" intiC (F)| for 0245F((J3). As shown 



in the figure, |-ffmt,c(F)| evolves upon cooling following 
the solid line M A fm{T) = M A fm(0)(1 - T/T N ) - 5 . 
Below T c =65 K, however, it is difficult to determine 
\Hi n t,c(F)\ because the spectral peaks are blurred, as 
shown in Fig. [7]Jc). The same situation is also true 
in 0245F(f|2); |i? int>c (F)| develops below T N =95 K as 
shown in Fig. [5Jd), but the F-NMR spectral peaks are 
blurred below T c ~ 75 K. Here, note that the Tjy values 
are expected to be the same in OP and IP because three 
dimensional AFM interactions determine T/y. In other 
words, we expect that Mafm(OP) and Mafm(IP) show 
the same T dependence, although F-NMR spectra probe 
the development of the staggered moment at OP below 
T N . 



3. F-NMR shift evidence for SC transition in AFM 
background 

Here, we deal with the SC properties in 0245F. The 
T dependence of 19 K c is displayed in Figs. IHtg)-|Hti) ■ 
The values of 19 K c are estimated from the gravity cen- 
ter of the 19 F-NMR spectra. In all samples, 19 K C is T- 
independent upon cooling down to T c , which is different 
from the Knight shift K S (T) of 63 Cu shown in Fig. H 
The reason that 19 K c is T-independent above T c is that 
the spin component in 19 K C is small owing to the small 
hyperfine coupling between 19 F nuclei and Cu-3d spins 
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as reported in the literature 64] . When decreasing T be- 
low T c , 19 K c markedly decreases due to the appearance 
of SC diamagnetism. Here, note that the reduction of 
19 K c , which is in association with the onset of high-T c 
SC, takes place under the background of the AFM order; 
these observations provide firm evidence for the uniform 
coexisting state of AFM and SC at a microscopic level. 

In 0245F(jj4), |f^ nt)C (F)| additionally increases below 
T c as shown in Fig. HJf). It is likely that the onset of 
high-T c SC decreases the size of (OP) due to the 

formation of coherent spin-singlet states over the sam- 
ple. Therefore, the additional increase in |-£/*m£,c(F)| may 
not be attributed to an increase of Mafm(OP) but an 
increase of the out-of-plane canting angle in the AFM- 
SC mixed state. In any case, this finding demonstrates 
an intimate coupling between the SC order parameter 
and Mafm [HI- The same behavior is expected for 
0245F(j}2) and 0245F(j}3) as well, although the F-NMR 
spectral peaks are blurred due to the spectral broadening 
in vortex states below T c . 



D. Phase diagram of AFM and SC in five-layered 
Ba 2 Ca 4 Cu50io(F,0)2 



phase, Tn listed in Table [TTI is plotted against p(IP). We 
assume that Tjy is determined by IP, considering the fact 
that Mafm (IP) is larger than M A fm (OP). The data 
point for Ijv=0 K at p ~ 0.126 corresponds to OP in 
0245F(til), which is paramagnetic even at T=1.5 K. 

This phase diagram demonstrates that in 0245F, the 
AFM metallic phase is robust up to p ~ 0.11 to 0.12, and 
that the three-dimensional long-range AFM order coex- 
ists with high-T c SC in an underdoped region. Such a co- 
existen ce p hase has been also reported in three-layered 
0223F Hi and four-layered 0234F [33|. Note that the 
coexistence of AFM and SC is not a phase separation 
between magnetic and paramagnetic phases; AFM and 
SC uniformly coexist in a Cu02 plane. We have dis- 
cussed the AFM ordering based on the zero-field NMR 
spectra shown in Fig. |6l where no paramagnetic signal is 
observed for Cu02 layers in AFM states. Here, we de- 
termine the critical hole density p c for the AFM order 
as p c ~ 0.11. The value of p c is expected to exist be- 
tween p(OP) ~ 0.106 in 0245F(j}2) and p(OP) ~ 0.126 
in 0245F((}1), because OP in 0245F(jjl) is paramagnetic, 
but OP in 0245F(fj2) is antiferromagnetic. 



IV. DISCUSSIONS 
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FIG. 9. (color online) AFM and SC phase diagram of five-layered 
0245F. Solid and open squares correspond to OP and IP, respec- 
tively. The values of p are determined from the room T values of 
Cu-NMR Knight shift in Fig. [4] and those of T N are from the 
T-dependence of [i?intc(P)| m Pig- El ( see t ex t for detail). The 
solid lines in the figure are guides to the eye. 

Figure [S] shows the phase diagram of AFM and SC 
in five-layered 0245F, which is derived from the present 
NMR study. For the SC phase, T c is plotted against 
p(OP). It has been reported that OP and IP have unique 
T c values, and that the higher one determines the bulk T c 
[52| . Here, it is expected in 0245F that the bulk T c listed 
in Table U corresponds to T c (OP). Both OP and IP are 
in underdoped regions, and p(OP) is larger than p(IP); 
therefore, T c (OP) is larger than T C (IP). For the AFM 



A. Layer number dependence of AFM-SC phase 
diagram in hole-doped cuprates 

Figure [TU] shows the phase diagram of hole-doped 
cuprates with the different stacking number n of Cu02 
layers. Figure [10] (a) is a schematic phase diagram for 
LSCO with ra=l; figs. CTbjlTTIe) show the phase dia- 
grams of Ba 2 Ca„_iCu n 2 „(F,0)2 (02(n-l)nF) with n=2 
[4l| , n=3 [Hj], n=4 [H, [65j], and n=5, respectively. The 
data in Fig.HOTe) correspond to those in Fig. HO and those 
in Figs.rTUTbVfTUrd) have been reported in previous NMR 
studies [23, [3a, |4l|. Here, note that the data of YBCO 
0, EH is cited as the AFM phase in Fig. [IWb). There 
are no data for the AFM phase of 0212F at present. 

Figures [TUTa'l- flW e') show the variation of the p c values 
for n-layered compounds, p c (n): p c (l}_~ 0.02 [lfj, l36j . 
vA?)_~ 0.055 0,113, Pc(3) ~ 0.075 (35|, p c (4) ~ 0.09 
[33|,|65|], andp c (5) ~ 0.11. This increase of p c {n) is qual- 
itatively explained as a result of the fact that the inter- 
layer magnetic coupling becomes stronger with increas- 
ing n. The mother compounds of high-T c cuprates are 
characterized by a large in-plane superexchange interac- 
tion Ji n ~ 1300 K between nearest-neighboring Cu-spins 
66 69]. In addition to ./,„. an Lnterlayer magnetic cou- 
pling along the c-axis plays a crucial role in stabilizing an 
AFM order since no long-range AFM order occurs in an 
isolated two-dimensional system at a finite temperature. 
The effective interlayer magnetic coupling of ri-layered 
cuprates is given as yj J c Jout{n), where J c is a magnetic 
coupling between OPs through CRL (charge reservoir 
layer), and J ut( n ) is that in a unit cell, as illustrated 
in Fig. UlT a). It is considered that J c is independent of 
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FIG. 11. (Color online) Schematic illustrations of magnetic cou- 
plings in n-layered cuprates. (a) Interlayer magnetic coupling along 
c-axis. J c is the magnetic coupling between OPs through CRL, 
which is independent of n; Jout(n) is the magnetic coupling in a 
unit cell, which increases with n. (b) In-plane superexchange in- 
teraction Ji n between nearest Cu spins in two-dimensional Cu02 
plane. The quantity Ji n is as large as 1300 K in undoped AFM- 
Mott insulators 66-69]. It is considered that Ji n does not de- 
pend on n; therefore, it is the effective interlayer magnetic coupling 
(y J c Jout(n)) that increases p c with increasing n (see text). 



FIG. 10. (color online) (a) Schematic phase dia- 

grams of LSCO. (b)-(e) AFM-SC phase diagrams of 
Ba 2 Ca n _iCu n O2n(F,O) 2 (02(n-l)nF). The data plot for (b) 
0212F, (c) 0223F, and (d) 0234F are cited from Ref. [3l 153. itl) . 
The data for (e) 0245F are the same with those in Fig. [9] The 
n dependence of the phase diagram reveals the variation of the 
p c (n) values: p c (l) ~ 0.02 [1(1 [3(|, Pc(2) ~ 0.055 0,0, p c (3) ~ 
0.075 [H, p c (4) ~ 0.09 [HHH, and p c (5) ~ 0.11. 



n, and that J ou t(n) increases with increasing n; therefore, 
the weakness of J ou t(n) suppresses the static long-range 
AFM orders in LSCO and YBCO at such small doping 
levels. 

While p c (n) increases with n, it seems to saturate at ~ 
0.14 to 0.16 even in the strong limit of interlayer magnetic 
coupling expected for infinite-layered compounds; Mafm 
at the ground state is extrapolated to zero at p ~ 0.14 to 
0.16, as discussed in the next section (see Fig. [T2|) . These 
results suggest that the uniform coexistence of AFM and 
SC is a universal phenomenon in underdoped regions, 
when the interlayer magnetic coupling is strong enough 
to stabilize an AFM ordering. A recent report seems to 
support this conclusion from a theoretical point of view 



B. Ground-state phase diagram of C11O2 plane 

Figure [T2ta) shows Mafm at T=1.5 K and the SC en- 
ergy gap Asc a s a function of p. The data include 0245F 
and theprevious reports for 0212F, 0223F, and 0234F 
[3314351 . l4lj |. The values of Mafm are estimated from 
zero-field NMR measurements, and those of Asc are as- 
sumed to be proportional to T c . We also plot Mafm 
for Y123 [HI in the figure instead of that for bi-layered 
0212F since the AFM phase of 0212F has not been re- 
ported. Note that, in contrast to Fig. [101 this figure 
focuses on the ground state of Cu02 planes, which gives 
us an opportunity to compare the experimental results 
with theoretical outcomes. 

The AFM phase in Fig. IT2T a) is characterized by the 
fact that Mafm exists up to the quantum critical hole 
density p qcp ~ 0.14 to 0.16, at which Mafm disappears 
even at the ground state. It is worth noting that irrespec- 
tive of n, the maximum of Asc (i-e., T c ) is around p ~ 
0.16, which is close to p qcp ~ 0.14. This implies an inti- 
mate relationship between the SC and AFM order param- 
eters. We also point out that at finite temperatures, no 
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FIG. 12. (color online) (a) Ground-state phase diagram of high-T c cuprates. The quantities of Mafm and Age are plotted as a function 
of p. The data for 0212F, 0223F, and 0234F are cited from Refs. [MUl El. and the data for Y123 are cited from Ref. 0. Here, it 
is assumed that Age is proportional to T c . The quantity Mafm exists up to the quantum critical hole density p qC p ~ 0.14 to 0.16, at 
which Mafm disappears even at the ground state. Note that the maximum of Age (i- e -. T c ) is around p ~ 0.16, which is close to p qC p- 
The inset shows the p-dependence of the supertransferred hyperfine fields B at Cu sites, which originates from the hybridization between 
Cu(3(i c 2_j / 2) and 0(2pa) orbits. The hybridization between Cu(3(l c 2_ H 2) and 0(2pa) orbits drastically increases as p increases above p 
~ 0.16 to 0.18, which would correspond to the disappearance of Mafm and the suppression of T c (sec text). The solid lines in the figure 
are guides to the eye. (b) AFM and SC phase diagram predicted by the t-J model [30j . M and 3> are the AFM and SC order parameters, 
respectively. 



static long range AFM order is observed when p c (n) < p, 
as shown in Fig. [TU1 although Mafm is expected to ex- 
ist at the ground state when p c (n) < p < p qcp - In the 
p range, AFM moments are fluctuating at finite temper- 
atures, which would induce various anomalies in under- 
doped regions. As a matter of fact, spin-glass phases or 
stripe phases have been reported in p c {ri) < P < Pqcp 
when n=l and 2 0, E3, S£ HI ; AFM moments fluctu- 
ating due to weak interlayer magnetic couplings may be 
frozen in association with disorders related to chemical 
substitutions, a buckling of Cu02 planes, and the onset 
of a charge order. 

As for the SC phase, Asc gradually decreases with 
decreasing p from the optimally-doping level, p ~ 0.16. 
The size of Asc (i-e., T c ) at p ~ 0.16 is largest when 
n=3, which is a universal trend in multilayered cuprates 
[7ll ufy . When we extrapolate the hole density p* at 
which A S c = 0, p* is ~ 0.05-0.06 for n—2; however, p* 
probably exists at below ~ 0.05 for rt=5. This variation 
of p* is clearly seen in Fig. [TU1 where it seems that p* 
moves to more underdoped regions as n increases. The 
reason for a possible reduction in p* is partly attributed 



to the flatness of the CuC>2 planes. With increasing n, it 
is expected that the flatness of CuC>2 planes is enhanced 
because the disorder effect related to heterovalent substi- 
tutions in CRL is relieved, especially at IP. Therefore, it 
would be possible that hole carriers slightly doped into 
Mott insulator can move even at low temperatures with- 
out Anderson localization when n=5, whereas they local- 
ize when n—1 and 2 [TH, |H, [73| . Further investigations 
on extremely underdoped regions for multilayered com- 
pounds remain as future works. 

We propose Fig. [12] as the ground state phase diagram 
of Cu02 planes. The ground state in underdoped regions 
is characterized by the uniformly mixed state of AFM and 
SC, which is observed even at finite temperatures when 
the interlayer magnetic coupling is strong enough to sta- 
bilize AFM orders. Another important point is that Asc 
(i.e., T c ) begins to decrease at p ~ 0.16 with increasing 
p. This is because a considerable change of electronic 
states occurs at p ~ 0.16. The inset of Fig. [T2l shows the 
p dependence of the supertransferred hyperfine field B 
at Cu sites in various hole-doped cuprates [4l|; B orig- 
inates from the hybridization between Cu(3c? a; 2_j / 2) and 
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0(2pa) orbits. As shown in the inset, the hybridization 
between Cu(3d a .2_ a 2) and 0(2pcr) orbits in overdoped re- 
gions are much larger than those in underdoped regions. 
This event corresponds to the suppression of the onsite 
Coulomb repulsion U or the in-plane superexchange in- 
teraction Ji n , which leads to the disappearance of Mafm 
and the suppression of T c . 

The ground state phase diagram presented here is qual- 
itatively consistent with theoretically predicted ones in 
terms of either the t-J model [IlM (2^, [H H, ,or 
the Hubbard model in a strong correlation regime [271 , 
[H[ . Figure [12jb) shows an AFM and SC phase diagram 
predicted by the t-J model with the condition that the 
second- nearest-neighbor hopping t' is zero [30j . Here, 
values for the vertical axis, M and $, in Fig. [TSTb) are 
an AFM and a SC order parameters, respectively. See 
Ref. [30J for details on t', M, and $. As shown in the 
figure, the AFM order parameter decreases with increas- 
ing p and vanishes at around p ~ 0.12 and coexists with 
the SC order parameter in underdoped regions. Our ex- 
perimental results demonstrate a good consistency with 
the t-J model; we conclude here that the large Ji n , which 
attracts electrons of opposite spins at neighboring sites, 
is the origin for high-T c SC. 



V. CONCLUSION 

Site-selective Cu- and F-NMR studies have unraveled 
the intrinsic phase diagram of AFM and SC in hole-doped 
cuprates Ba 2 Ca„_iCu„0 2 „(F,0) 2 (02(n-l)nF). The ob- 
tained results are as follows: 

(i) The AFM metallic state is robust up to the critical 
hole density p c for AFM orders, and it uniformly 
coexists with SC. 

(ii) The p c for n-layered compounds, p c (n), increases 



with increasing n due to the growth of the interlayer 
magnetic coupling. The p c (n) values are deduced 
asp c (l) ~ 0.02 0,l23,p c {2)~ 0.055 0,03, Pc(3) 
- 0.075 [Hj], p c {4) ~ 0.09(33, and Pc(5) ~ 0.11. 

(iii) The maximum of A.$c (he., T c ) takes place around 
p ~ 0.16, which is just outside the quantum critical 
hole density p qC p ~ 0.14 to 0.16. 

(iv) The ground-state phase diagram of AFM and SC 
is in good agreement with that theoretically pre- 
dicted by the t-J model [TIM [H, H H, \M or 
by the Hubbard model in the strong correlation 
regime [H |Hj]. 

We conclude that those results are accounted for by the 
Mott physics based on the t-J model. In the physics be- 
hind high-T c phenomena, there is a very strong Coulomb 
repulsive interaction U (>6 eV), which prohibits the dou- 
ble occupancy of an up-spin electron and a down-spin 
electron at the same site. When noting that the strength 
of U is almost unchanged by hole doping, it is considered 
that the large Ji n attracts electrons of opposite spins at 
neighboring sites [13, [H, [29| . The qualitative consistency 
in the ground-state phase diagram between our experi- 
mental results and theoretical ones support the t-J model 
as a mechanism of high-T c SC. 
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